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We study a random magnetic resonance force microscopy (MRFM) signal 
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\ caused by the thermal vibrations of high frequency cantilever modes in the 



oscillating cantilever-driven adiabatic reversals (OSCAR) technique. We show 
that the regular MRFM signal with a characteristic decay time, r m , is followed 
by a non-dissipative random signal with a characteristic time r r . We present 
the estimates for the values of r m and r r . We argue that this random MRFM 
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signal can be used for spin detection. It has a "signature" of a sharp peak in 
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. its Fourier spectrum. 
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& ! I. INTRODUCTION 
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Magnetic resonance force microscopy (MRFM) is a powerful combination of magnetic 

O ' 

resonance and atomic force microscopy which promises an atomic resolution for the magnetic 
•^H ■ resonance studies [1,2]. When the number of detecting spins is small enough the fluctuating 
magnetic moment of a spin system becomes comparable with the maximum possible mag- 
netic moment of the system. In this case, one expects to observe a random MRFM signal 
which follows a regular MRFM signal. This paper is devoted to the first theoretical analysis 
of this random MRFM signal. We consider the random signal caused by the high frequency 
thermal cantilever vibrations. 

Theoretical analysis for the initial (regular) MRFM signal in the OSCAR technique was 
performed in [3,4]. In [3] the Bloch type equations were derived to describe the motion of the 
local magnetization in the "resonance slice" of a sample. In these equations the relaxation 
terms are associated with the fluctuating magnetic field caused by the thermal vibrations of 
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the cantilever. The approach used in [3] implies no spatial correlations for the fluctuating 
magnetic field. As a result, the local magnetization in the resonant slice monotonically 
decays. In [4] we have developed a computational model for simulations the MRFM signal 
in the OSCAR technique. Numerical simulations allow us to study both regular and random 
OSCAR MRFM signals. In this paper we use computer simulations to analyze the random 
MRFM signal. 

In Section 2, we briefly explain the OSCAR technique, in particular, the sign of the 
frequency shift of the cantilever vibrations. In Section 3, we describe our numerical model. 
In Section 4, we present the qualitative analysis of the spin motion under the action of the 
fluctuating magnetic field caused by the high frequency cantilever vibrations. In Section 
5, we describe the results of our numerical simulations of the regular and random MRFM 
signals. 

II. SYSTEM UNDER CONSIDERATION 

We consider the "parallel setup" used in the experiment [2]. (See Fig. 1.) We study the 
oscillating cantilever-driven adiabatic reversals (OSCAR) technique which is currently the 
most promising MRFM technique for spin-based atomic scale resolution [2,5]. In OSCAR 
a feed-back loop is used to maintain a constant amplitude of the cantilever vibrations. 
The cantilever vibrations, in combination with an external rf magnetic field, cause cyclic 
adiabatic reversals of the magnetic moment in the tiny "resonance slice" of the sample. 
In turn, the magnetic reversals generate an oscillating magnetic force on the cantilever tip 
which is approximately proportional to the cantilever amplitude. This force changes the net 
restoring force and, consequently, changes the frequency of the cantilever vibrations which 
is to be detected. 
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FIG. 1. The "parallel setup" used in the experiment [2]. B ext and Bi(t) are the permanent and 
rotating rf magnetic fields; rh is the magnetic moment of the ferromagnetic particle; fa is the fc-th 
magnetic moment in the resonant slice. 

The sign of the frequency shift can be easily understood from Fig. 1. Consider, for 
simplicity, only the magnetic moment fa at the center of the resonant slice. When the 
z coordinate of the ferromagnetic particle z c is positive, the distance between the particle 
and the center of the resonant slice is greater than the equilibrium distance. Consequently, 
the z component of the magnetic field at the center of the resonant slice is less than the 
same quantity in the equilibrium position z c = 0. Then, in the rotating reference frame, 
the z component of the effective field at the center of the resonant slice is negative. If the 
magnetic moment fa points in the direction of the effective field (the ground state), then the 
z component of the magnetic moment ji z is also negative, unlike the magnetic moment of the 
cantilever tip. Because opposite magnetic moments repeal each other, the magnetic force 
on the ferromagnetic particle points in the positive z direction, opposite to the cantilever 
spring force. As a result, the cantilever frequency decreases. Based on the same arguments, 
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if the magnetic moment in the resonant slice is opposite to the direction of the effective field, 
its magnetic force causes the cantilever frequency to increase. 



III. THE MODEL FOR NUMERICAL SIMULATIONS 

We describe the vibrations of the cantilever tip or, more precisely, the center of the 
ferromagnetic sphere, as the motion of a harmonic oscillator with a fundamental frequency 
oj c and a quality factor Q: 

z c + z c + z c /Q = /(r). (1) 

In this equation, we use the dimensionless coordinate z c of the center of the ferromagnetic 
particle (in terms of the amplitude of the cantilever vibrations z ), and the dimensionless 
time r = u c t, /(r) is the dimensionless magnetic force produced by the spins of the resonant 
slice on the ferromagnetic particle: 

N 

= ^2VkVkz, (2) 
k=l 

Ho 3m/j,z k (5z 2 k - 3rl) _ 

% = =7 > r >< = ( x k,yk,z k ). 

k c z r k 

Here N is the number of magnetic moments in the resonant slice, fl k is the fc-th magnetic 
moment (in terms of its magnitude fi k ), k c is the effective spring constant of the cantilever, m 
is the magnetic moment of the ferromagnetic particle, f k = f k (z c ) is the vector connecting 
the center of the vibrating ferromagnetic sphere with the position of the k-th magnetic 
moment. 

We assume the number of spins N in the resonant slice is large enough, so that we can 
ignore quantum effects and consider classical equations of motion for the magnetic moments 
jl k , in the rotating frame: 

Afcrr = — Afc/ifc^, (3) 

fi ky = A k fj, kx — e/i kz , 
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Here 

A fc = (<yB ext -u ) uj c +^- 3 fc fc , 4 

47r cu c Zo r£ 

7 is the magnitude of the spin gyromagnetic ratio, u is the frequency of the rf field. 

To simulate the thermal noise caused by the high frequency cantilever modes, we replace 
the vector r k (z c ) in the expression for A k by the vector r k [z c + Sz c ), where 

5z c = J2 2a nCos(tt n T + tf n ). (5) 

n 

In Eq. (5), a n and Q n are the dimensionless thermal amplitude and the frequency of the 
n-th cantilever mode, and \l/ n is a random phase. The factor 2 in Eq. (5) appears because 
the amplitude of the cantilever tip for any mode n is twice the amplitude of the mode [6]. 
The thermal amplitude of a high frequency mode can be estimated from the equipartition 
theorem 

a n = (z Q n )~ 1 (k B T/2k c ) 1/2 , (6) 

where T is the cantilever temperature. In Eq. (6) we took into account the expression 
for the cantilever fundamental mode frequency ul = 4k c /m c , where m c is the mass of the 
cantilever. (See, for example, [7].) 



IV. QUALITATIVE ANALYSIS AND ESTIMATES 

Consider the motion of a single classical magnetic moment jx k in the center of the resonant 
slice, under the action of the fluctuating magnetic field produced by the ferromagnetic 
particle. Assume that the magnetic moment moves adiabatically, together with the effective 
field in the semi-plane (+z)-(+x)-(-z), in the rotating frame. When the polar angle of the 
vector /2fc is not small, the weak random fluctuating field produced by the ferromagnetic 
particle has a component perpendicular to /4 which causes a deviation from the effective 
field. The resonance frequency u) e in the rotating frame is 

„ e = K 2 + ^) 1/2 , (7) 
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where u z = oj z (z c ) is associated with the z component of the effective field which is pro- 
portional to the ferromagnetic particle coordinate z c , and uj r = r yB 1 is the Rabi frequency 
associated with the rf field. The maximum value of u z , u zm , is much greater than ujr. In 
the process of adiabatic reversals the resonance frequency of the magnetic moment changes 
from uj zm (near the (+z) axis) to u R (near the transverse plane), and back to uj zm (near the 
(-z) axis). If a characteristic frequency of the fluctuating field falls in the region {u Rl oj zm ) 
the fluctuating field causes a noticeable deviation of the magnetic moment from the effec- 
tive field. We consider the fluctuating field generated by high frequency cantilever modes 
which cover the interval (u>R,iv zm ). The thermal amplitudes of these modes are inversely 
proportional to their frequency. (See Eq. (6).) As an example, Fig. 2 demonstrates the 
thermal amplitudes of the cantilever modes in the interval (cur, 5u r ) for the data presented 
in experiment [2] . One can see that among the high frequency modes which can be resonant 
to the magnetic moment, the greatest thermal amplitudes have the modes in the narrow 
region, approximately (u R ,2u R ). The geometrical factor also favors these modes because 
the magnetic moment jlk has the resonant frequency in the region (cur, 2u r ) when it is very 
close to the transverse plane where the fluctuating field is perpendicular to the magnetic mo- 
ment /2fc. To estimate the action of the fluctuating field, we will consider only the modes in 
the narrow region (u>r, 2ur). We will assume that in this region of the resonant frequencies 
the magnetic moment experiences the fluctuating magnetic field produced by the thermal 
vibrations of the cantilever. To simplify our estimate, we assume that these vibrations have 
the constant amplitude a-r and the random phase. We estimate the amplitude ax using Eq. 
(6) and substituting Q n — > uj r /uj c : 

a T = (co c /u JR )(k B T/2k c ) 1 / 2 . (8) 

The corresponding thermal amplitude of the ferromagnetic particle on the cantilever tip is 
2clt- When the magnetic moment passes the region of the resonant frequencies (u R ,2u R ), 
it experiences a deviation from the effective field. The amplitude of the fluctuating field 
is Gclt where G is the magnetic field gradient, and factor 2 disappears as we consider the 
rotating component of the fluctuating field. 
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The characteristic time between the "phase jumps" of the fluctuating field (the correla- 
tion time r ) can be estimated as the reciprocal of the frequency region r pa 1/ujr. Assuming 
that the angular deviation from the effective field is a diffusion process, we obtain 

A©f = DAti, (9) 

where AO 2 is the characteristic deviation during a single reversal, D is the diffusion coef- 
ficient, and Aii is the time of passing the region (u>r, 2ujr) during one reversal. The value 
of D can be estimated as A0q/t o , where A6 is the characteristic deviation between two 
consecutive phase jumps, A0 O = ^GcltTo. The time interval Aii can be easily estimated 
from the equation u e = 2ujr. Substituting into Eq. (7) 

u) z = r yGz cosu c t, (10) 

we obtain 

Aii = 2y/3(*yGzo)- 1 (u R /u c ). (11) 

(Note that the number of phase jumps during one reversal can be estimated to be Aii/r = 
ujRAti. For example, for experiments [2], where uj c /2tt = 21.4 kHz, G = 1.4 x 10 5 T/m, 
z = 28 nm, ujr/2tc ranges between 4.17 MHz and 8.34 MHz, the number of jumps ranges 
between 25 and 100.) Combining Eqs. (9) and (11) we find 

A0f = (2V3)(-fGa 2 T /uj c zo). (12) 

Now we can estimate how many reversals are needed to decrease twice the component of 
the magnetic moment along the effective magnetic field. To get the total deviation A© 2 pa 1 
we need 1/A0 2 reversals. The corresponding time r m is 

r m = 7r/(cu c A0f) = (7r/2v / 3)(^o/7G'a|). (13) 

Substituting in Eq. (13) the expression for a T in (8), and omitting the numerical factor 
ir/y/3, we obtain finally the required estimate 

r m ~ {k c z G hGk B T){u R lu c ) 2 . (14) 
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We shall discuss now this formula. The characteristic time of the deviation from the 
effective field increases with increasing Rabi frequency cur, but decreases with increasing 
temperature. These dependences reflect the obvious properties of the thermal noise of the 
cantilever. The dependence r m on the amplitude of the cantilever vibrations zo is associated 
with the time Ati of passing the resonant region {oj r ,2uj r )\ the greater z is, the smaller 
time Ati, thus the greater number of reversals is needed to provide a significant deviation 
of /4 from the effective field. The dependence r m on the magnetic field gradient G appears 
as the result of two competing factors. On the one hand, when the gradient G increases, 
the fluctuating magnetic field also increases. On the other hand, the time of passing the 
resonant region Ati decreases when the gradient G increases. 

Now consider the motion of the magnetic moment in the "tilted system" connected to the 
effective field. In the "titled system" a single classical magnetic moment moves randomly 
between the direction of the effective field and the opposite direction. If the characteristic 
time of the deviation from the effective field is r m , then after time ~ 2r m the MRFM signal 
changes from its maximum value Au m to zero. After this, the MRFM signal may change 
randomly with the characteristic time r r ~ 8r m , sometimes approaching its maximum value 
Au m . (Correspondingly, a single quantum spin 1/2 is expected to jump randomly between 
two directions of the effective magnetic field with the characteristic time r r ~ 8r m . In the 
quantum case, the MRFM signal will take only two values ±Auj m .) For a large number 
of spins in the resonant slice, the quantum and classical systems have a similar behavior: 
one may expect first a "regular decay" of the MRFM signal with the characteristic time 
r m , and then the random MRFM signal of a smaller amplitude with the characteristic time 
r r ~ 8r m . The amplitude of the random signal is expected to be smaller than the maximum 
value of the regular signal Au m because spins in the resonant slice change their directions 
at different times. So, they may experience different fluctuating fields and have different 
components along the effective field. 

Finally, if the number of spins in the resonant slice is very large, one can choose a small 
(but still macroscopic with many spins) volume of averaging AV, where all spins experience 
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the same effective field (the external field and the field produced by the ferromagnetic 
particle in the rotating frame). After this, one can introduce a continuous magnetization in 
the resonant slice 

where the sum is taken over the spins in the volume AV^ near the point r. Note that this 
approach assumes that the effective field inside the volume of averaging is approximately 
uniform. So, the magnitude M is constant. To describe the change of the magnitude 
M due to the microscopic relaxation mechanisms (e.g. the spin-phonon interaction) one 
adds the relaxation terms in the equations of motion, e.g. Bloch type relaxation terms. 
In our case, the fluctuating field is produced by the ferromagnetic particle. The scale of 
the spatial inhomogeneity for the fluctuating field is the same as for the regular magnetic 
field produced by the ferromagnetic particle. Thus, we expect that the local magnetization 
magnitude M{r) will conserve until the microscopic relaxation mechanisms come into play. 
As a result, we may expect that the random MRFM signal will be observed even for a 
continuous magnetization. 

V. NUMERICAL SIMULATIONS 

In this section, we present the results of computer simulations of the random MRFM 
signal in the OSCAR technique. We solved Eqs. (1) and (3) with thermal noise (5). The 
magnetic moments were uniformly distributed in the resonant slice. The upper and the 
lower boundaries of the resonant slice have been determined from the equations = 
at z c = ±1. Each time the cantilever passes the upper point, the computer automatically 
increases the value of z c to 1. This procedure corresponds to the action of a feed-back loop 
in OSCAR which maintains a constant amplitude of the cantilever oscillations. The main 
parameters of the problem were taken from the experiment [2] 

B ext = 140mT, k c = 0.014N/m, oo c /2ir = 21.4kHz, z = 28nm, Q = 2 x 10 4 , 

m = 1.5 x 1(T 12 J/T, M = 0.89A/m, d 1 = d/2 = 700nm, d 2 = 875nm, G = 1.4 x 10 5 T/m, 
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where d is the diameter of the ferromagnetic particle. The frequencies Vt n have been calcu- 
lated for a silicon rectangular cantilever. (See Fig. 2.) 
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FIG. 2. The thermal amplitudes of the high frequency cantilever modes for a silicon cantilever 
(190^m x 3/xm x 850nm) used in [2]. The left-most bar corresponds to the seventh harmonic. The 
Rabi frequency oor/2tt = 4.17MHZ. The cantilever temperature is 80K. 

We calculate the relative shift of the cantilever period AT/AT caused by the spins of 
the sample (where AT is the initial shift of the cantilever period). Initially all magnetic 
moments point approximately in the direction of the effective magnetic field. So, the initial 
shift AT is positive. The phases \I/ n in Eq. (5) were changed randomly in the interval 
(0, 27r). The characteristic time interval between the phase jumps was taken as 10 Rabi 
periods 10 x (2tt/ur). The number of high frequency modes taken into consideration was 
22. 

As an example, Fig. 3 demonstrates the connection between the regular and random 
MRFM signals, for a system of 50 magnetic moments uniformly distributed in the resonant 
slice. 
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FIG. 3. The random MRFM signal. The insert shows the regular MRFM signal. The tempera- 
ture is 80K. The Rabi frequency ujr/2tt = 4.17MHz. The cantilever amplitude zq = 28nm. 




FIG. 4. The Fourier spectrum of the random MRFM signal shown in Fig. 3. 
The decay time of the regular signal is r m = 145 ms. (The value of r m estimated 
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using the formula (14) is 550 ms. Thus, our rough estimate provides the same order of 
magnitude as the computer simulations.) The MRFM random signal is about 30% of the 
maximum regular signal, and the characteristic time of the random signal r r = 895ms is not 
far from 8r m . Fig. 4 shows the Fourier spectrum of the random signal (in arbitrary units): 
AT(t)/AT = J2 A m cos(is m t + $ m ). One can see a sharp peak near the value l/r r . In fact, 
we find the value r r from the position of the peak in the Fourier spectrum. This position is 
a "signature" of the random MRFM signal. 

Fig. 5 shows the standard deviation of the random signal = \J £ — £) 2 (where £(t) is 
the random function £(t) = AT(t)/ AT ) as a function of a number of magnetic moments 
N distributed in the resonant slice (for the same value of the average magnetization M). 
One can see that with increasing N the standard deviation approaches the value 0.1. 
This indicates that the random MRFM signal survives the transition to the continuous 
magnetization. Note that the position of the peak in the spectrum of the random signal (see 
Fig. 4) does not depend on N. 




FIG. 5. The standard deviation of the random MRFM signal as a function of the number N 
of magnetic moments (cells) in the resonant slice (at a fixed value of the average magnetization 
M = 0.89A/m). All parameters are the same as in Fig. 3. 
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FIG. 6. The component [i\\(t) of an arbitrarily chosen magnetic moment along the effective 
magnetic field in the rotating frame. All parameters are the same as in Fig. 3. 

As an illustration to the spin dynamics, Fig. 6 shows the random change of the magnetic 
moment component along the effective field. (The magnetic moment was arbitrary chosen 
in the resonant slice.) One can see that the magnetic moment randomly moves between the 
direction of the effective field and the opposite direction. The characteristic time of this 
motion r r is the same as for the MRFM signal shown in Fig. 3. 

As an extreme example, we have considered the case when the cantilever amplitude is 
only 60 pm. In this case, the initial magnetic moment of the resonant slice for our value 
of the magnetization M = 0.89 A/m, is approximately one Bohr magneton. Figs. 7 and 
8 show the regular and random signals for the number of magnetic moments (cells in the 
resonant slice) N=200 and N=1000. One can see that the relaxation time of the regular 
signal r m 5.5 ms is the same for both values of N. 



13 




0,0-1 . 1 ■ 1 . 1 . 1 . 1 

2 4 6 8 10 

t, ms 



FIG. 7. A regular MRFM signal for an extremely small amplitude of the cantilever vibrations, 
zq = 60pm. The temperature of the cantilever is 0.57K. The Rabi frequency oor/2tt = 8.34MHz. 
Curve 1: N=200; curve 2: N=1000. 
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The value of the standard deviation for the random signal drops from 0.105 to 0.052, 
when N increases from N=200 to N=1000. The search for a "plateau" in the dependence, 
a^(N), requires a greater number of the magnetic moments (cells) and, correspondingly, 
much greater computational time. The spectrum of the random MRFM signal for both 
values of N has its maximum at a frequency of approximately 17.5 Hz, which corresponds to 
the characteristic time r r = 57 ms. This is not far from the estimated value r r ~ 8r m = 44 
ms. Note, that for the extremely small amplitude of the cantilever vibrations considered here, 
instead of adiabatic reversals we have small adiabatic oscillations of the effective magnetic 
field near the x-axis in the rotating frame. 

CONCLUSION 

In this paper we study the random non-dissipative MRFM signal caused by the high 
frequency thermal vibrations of the cantilever. For a single spin detection, this random 
signal is expected to have the same amplitude as the initial regular signal. We show that for 
a relatively large number of spins in the resonance slice the random non-dissipative signal 
occurs after the relatively short regular signal. Our computer simulations indicate that 
even for a very large number of spins, when one can use the approximation of a continuous 
magnetization, a long lasting random signal can be observed, until the intrinsic mechanisms 
of spin relaxation destroy the local magnetization in the resonant slice. The random non- 
dissipative MRFM signal can be used for the spin detection, including a single spin detection. 
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